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ABSTRACT 

We have used the ATNF Mopra antenna and the SEST antenna to search in the 
directions of several class II methanol maser sources for emission from six methanol 
transitions in the frequency range 85-115 GHz. The transitions were selected from 
excitation studies as potential maser candidates. Methanol emission at one or more 
frequencies was detected from five of the maser sources, as well as from Orion KL. 
Although the lines are weak, we find evidence of maser origin for three new lines in 
G345.01+1.79, and possibly one new line in G9.62+0.20. 

The observations, together with published maser observations at other frequencies, 
are compared with methanol maser modelling for G345.01+1.79 and NGC 6334F. We 
find that the majority of observations in both sources are consistent with a warm 
dust (175 K) pumping model at hydrogen density ~ 10 6 cm~ 3 and methanol column 
density ~5x 10 17 cm~ 2 . The substantial differences between the maser spectra in the 
two sources can be attributed to the geometry of the maser region. 

Key words: Hn regions — ISM: molecules — masers — radio lines : ISM 



1 INTRODUCTION 

Class II methanol masers provide a unique high resolution 
tracer of physical conditions at the sites of massive star for- 
mation. The two strongest class II methanol maser tran- 
sitions are the 5i — 6o A + transition at 6.6 GHz and the 
2o — 3_i E transition at 12.1 GHz. Single dish and inter- 
ferometric studies of these transitions (e.g. Menten 1991, 
Caswell et al. 1995a, b, Norris et al. 1993) have revealed the 
widespread distribution of the masers, the forms of maser 
spot distributions, and their velocity fields. However, data 
on only two transitions are not sufficient to determine phys- 
ical parameters in the sources through modelling because 
there are so many model parameters (see, e.g., Sobolev, 
Cragg & Godfrey 1997a). 

In addition to the very bright 6.6- and 12.1-GHz 
methanol masers, moderately bright class II masers have 
been identified at 19.9 GHz (Wilson et al. 1985), 23.1 GHz 
(Wilson et al. 1984), 28.9 GHz (Wilson et al. 1993), 37.7, 



38.2 and 38.4 GHz (Haschick, Baan & Menten 1989), 
107.0 GHz (Val'tts et al. 1995), 156.6 GHz, and in a se- 
ries of eight lines near 157 GHz (Slysh, Kalenskii & Val'tts 
1995). More recently a new weaker class II maser line has 
been identified at 108.8 GHz (Val'tts et al. 1999). 

A model featuring pumping by IR radiation to the sec- 
ond and first torsionally excited states of methanol, together 
with amplification of background H II region continuum ra- 
diation, can account for the extreme brightness of the 6.6- 
and 12.1-GHz masers (Sobolev & Deguchi 1994a, Sobolev 
et al. 1997a). Furthermore, it predicts maser action span- 
ning a wide range of brightness temperatures in ~ 100 other 
methanol lines (Sobolev, Cragg & Godfrey 1997b, hereafter 
SCG97). This list of maser candidates includes all 18 estab- 
lished class II masers and also the more recently detected 
108.8-GHz weak maser. 

Here, we extend the search for weaker class II methanol 
maser lines, based on the list of maser candidate transitions 
from SCG97. Sections 2 and 3 report a search for the five 
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Table 1. Methanol laboratory frequencies from Tsunekawa et al. 
(1995). Figures in brackets represent uncertainty in kHz. The last 
column is the velocity resolution (after Hanning smoothing) of the 
astronomical observations reported here. 



Vt 


Transition 


Frequency 


Resolution 






(MHz) 


(kms -1 ) 





6-2-7-iE 


85568.084(30) 


0.28 





7 2 - 6 3 A~ 


86615.578(30) 


0.43 





7 2 - 63 A+ 


86902.956(50) 


0.43 


1 


lo - 2i E 


93196.682(30) 


0.40 





83 - 9 2 E 


94541.778(50) 


0.40 





7 2 - 81 A+ 


111289.515(50) 


0.34 



remaining transitions in the 85-115 GHz band not previously 
observed in class II maser sources. We also searched for a 
line of the first torsionally excited state of methanol (vt = l), 
predicted to have class II maser character, but too weak 
to be included by SCG97. The quantum numbers and rest 
frequencies of the six lines are given in Table [l]. 

Since the relative brightness of different methanol 
masers is sensitive to the model parameters (Section 4.1), 
multi-transition studies provide a probe of the physical con- 
ditions in star-forming regions if the masers are co-spatial 
(Section 4.2). We compare the model with all available 
methanol maser observations for G345.01+1.79 (Section 4.3) 
and NGC 6334F (Section 4.4). 



2 OBSERVATIONS 

The 85.5-GHz observations were made at SEST0 1996 
March using standard SEST beam switching and calibration 
(tittp://www. Is. eso.org/lasilla/Telescopes/SEST). The high 
resolution acousto-optical spectrometer provided a 2048- 
channel spectrum over an 86-MHz bandwidth centred near 
85568 MHz. The antenna half-power beam width (HPBW) 
at 85.5-GHz is 58", the conversion efficiency for a point 
source is 25.2 Jy/K, and the pointing accuracy is better than 
5". The system temperature was about 130 K, and signal in- 
tegration periods were 12 min (exceptions were G9. 62+0. 20 
- 32 min and G12.9T0.26 - 8 min). The data were processed 
using 'CLASS' from the GILDAS software package. 

The remaining observations were made with the Mopra 
Antenna, operated by the Australia Telescope National Fa- 
cility, ATNf|} Observations at 86.6, 86.9 and 93.1 GHz were 
made in 1996 November, and at 94.5 and 111.2 GHz in 1998 
May. The illuminated diameter of the antenna was 15.4 m, 
with HPBW 37" at 111 GHz increasing to 48" at 86 GHz. 

Spectra from a digital correlator comprised two bands, 
each of 1024 channels spanning 64 MHz. One band was 
tuned to the methanol transition, giving a channel sepa- 
ration of approximately 0.2 kms -1 (Table[j|). Generally the 
second band was reserved for 86.2-GHz SiO maser point- 
ing calibration, but the 86.6-GHz methanol observation was 



* SEST=Swedish-ESO Submillimetre Telescope, La Silla, Chile 
1 The Australia Telescope is funded by the Commonwealth of 
Australia for operation as a National Facility managed by CSIRO. 



made simultaneously with the 86.9-GHz observation by re- 
tuning this second band. On-source observations were inter- 
leaved with reference spectra (offset by 30' or 60' in decli- 
nation) using five-minute integration times. 

The 1996 November data at 86.6, 86.9 and 93.1 GHz 
were not systematically corrected for atmospheric attenua- 
tion; rather, the recorded antenna temperatures have been 
scaled by a factor of 2, indicated from comparison between 
our 86-GHz Orion KL spectra and those observed by Jo- 
hansson et al. (1984) at Onsala, and by Turner (1991) at 
Kitt Peak, while taking into account the different beam 
sizes. The intensity calibration is uncertain to a factor of 2, 
but adequate for a detection experiment, and for order-of- 
magnitude comparisons with modelling. We applied a con- 
version factor of 50 Jy/K to the corrected antenna temper- 
atures. 

During the 1998 May observations at 94.5 and 
111.2 GHz, additional measurements of receiver total power 
were made with an ambient temperature load (assumed tem- 
perature 290 K) regularly inserted. These alloow correction 
for atmospheric absorption using the method of Kutner & 
Ulich (1981). The uncertainty of the absolute flux density 
scale of these observations is 20%. 

The data were processed using the SPC data reduction 
package at the ATNF. After averaging individual spectra 
and subtracting baselines, Hanning smoothing was applied 
and Gaussian line profiles fitted. Velocity tracking was not 
available at Mopra, but the LSR velocity shift over the time- 
frame of our observations was generally less than one chan- 
nel, so that the spectra are not significantly smeared. A fre- 
quency setting problem at 93.1 GHz led us to adjust the 
velocity scale by 1.61 kms -1 following comparison of the 
sharply defined spectrum of N 2 H + in L134N with previous 
observations (Womack, Ziurys & Wyckoff 1992). 

Class II maser sources selected from the surveys of 
Caswell et al. (1995a, b) are listed in Table ^| and were ob- 
served at the various methanol frequencies as listed in Ta- 
ble [l]. As a system test at Mopra we reobserved previously 
detected thermal methanol emission from Orion KL. 



3 RESULTS 

The present observations detected weak high frequency 
methanol emission from five of the class II methanol maser 
sources, as well as Orion KL. The parameters of peak flux 
density, mean velocity, and width to half intensity, as derived 
from Gaussian fits to the observed methanol lines, are given 
in Table []. Upper limits where no line was detected are sum- 
marised in Table ^| The spectra of detected methanol lines 
are presented in Figs. 1-4. 

It is not always easy to distinguish between weak maser 
emission from a very compact source and thermal emission 
from an area that is more extended, but nevertheless smaller 
than a single dish telescope beam. Many sources of class II 
methanol maser emission also exhibit thermal methanol fea- 
tures at millimetre wavelengths (e.g. Caswell et al. 2000). 
In identifying possible new masers we have been guided by 
the width of the line and the velocity range exhibited by 
previous maser and thermal observations. In the sources ex- 
amined all firmly established maser features have widths of 
the order 0.5 kms -1 , while the thermal features are con- 
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Table 2. Upper limits in Jy for methanol line search (3 X rms); D denotes detection, with details in Table |3j. Source coordinates are 
from Caswell et al. (1995a). 
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G336.02-0.83 


16 


31 


24.8 


-48 


10 


37 




<2.7 


<4.6 








G339.88-l.26 


16 


48 


24.76 


-46 


03 


33.9 




<1.7 


<2.1 


<5.3 


<2.5 




G345.00-0.22 


17 


01 


38.5 


-41 


24 


59 




<2.9 
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G345.01+1.79 


16 


53 
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09 


46.0 
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NGC 6334F 


17 


17 


32.3 
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14 
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D 
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G351. 77-0.54 


17 


23 


20.67 


-36 


06 


45.4 
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D 








G9.62+0.20 


18 


03 


15.98 
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52.9 


D 


<2.3 
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G10.47+0.03 
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05 
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G12.68-0.18 
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10 
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<5.2 








G12.91-0.26 
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11 


43.8 


-17 


53 


04 


<3.8 












G23.01-0.41 
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31 


55.6 
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03 


09 


<3.8 












G31.28+0.06 
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45 


37.2 


-01 


30 


00 


<3.8 












G35.20-l.74 


18 


59 


13.1 


+01 09 07 


<3.8 


<2.3 


<2.8 


<4.9 


<3.6 





siderably broader in the range 3-10 kms" 1 . Only fea- 
tures narrower than 2 kms" 1 appearing at the character- 
istic methanol maser velocity of the source are considered 
here as probable new masers (see discussion of sources in 
the following subsections). 

3.1 Orion KL 

Strong methanol emission was detected at all four frequen- 
cies observed in Orion KL (Fig. |l|) at flux density levels 
exceeding any of the maser sources. The lines have all been 
previously observed (Lovas et al. 1976, 1982, Hollis et al. 
1983, Johansson et al. 1984, Turner 1991) and are attributed 
to thermal emission at rotational temperature 100-200 K. 

The general shape of the Orion KL spectra in the four 
transitions is similar to the 107- and 108-GHz methanol 
spectra (Val'tts et al. 1999) but weaker. Val'tts et al. inter- 
pret their 107- and 108-GHz spectra as thermal, but with 
broad and narrow components, both at a velocity of ap- 
proximately 7.5 kms. Our spectra show less evidence for any 
narrow component. 

3.2 G345.01 + 1.79 

Narrow methanol lines were observed at 85.5, 86.6 and 
86.9 GHz near velocity -22 kms" 1 in G345.01+1.79 (Fig.|). 
These coincide with strong maser peaks at -21.8 kms" 1 in 
the 6.6- and 12.1-GHz spectra (but note an even stronger 
6.6-GHz peak at -17.9 kms" 1 ). Previous methanol obser- 
vations of this source at higher frequencies display narrow 
maser emission near -22 kms -1 accompanied by (but clearly 
distinguishable from) broad thermal emission centred near 
-14 kms" 1 . Such behaviour is seen at 156.6 and 107.0 GHz 
(Caswell et al. 2000), at 108.8 GHz (Val'tts et al. 1999) and 
in the 157-GHz series (Slysh et al. 1995), and is also il- 
lustrated in Fig. for the 108.8-GHz emission. The thermal 



Table 3. Summary of methanol line detections. Figures in brack- 
ets are uncertainties from Gaussian line fitting. 



Source S V LS r AV 

(Jy) (kms" 1 ) (kms" 1 ) 



85.5 GHz 








G345.01+1.79 


10.0(0.8) 


-22.1(0.1) 


0.8(0.1) 


NGC 6334F 


7.7(1.0) 


-8.3(0.1) 


5.4(0.1) 


G9.62+0.20 


0.9(0.1) 


3.8(0.3) 


5.7(0.6) 


// 


1.2(0.4) 


-0.9(0.1) 


0.5(0.1) 


G10.47+0.03 


3.3(1.0) 


66.2(0.2) 


9.0(0.4) 


86.6 GHz 








Orion KL 


12.3(0.7) 


8.0(0.3) 


4.2(0.3) 


G345.01+1.79 


2.8(0.3) 


-21.9(0.3) 


1.5(0.2) 


NGC 6334F 


2.6(0.2) 


-8.1(0.4) 


5.5(0.4) 


G351. 77-0.54 


1.5(0.2) 


-3.1(0.6) 


9.8(1.0) 


86.9 GHz 








Orion KL 


14.2(0.3) 


8.1(0.3) 


4.4(0.1) 


G345.01+1.79 


4.1(0.4) 


-21.7(0.3) 


1.0(0.2) 


NGC 6334F 


3.6(0.2) 


-8.0(0.3) 


4.6(0.3) 


G351. 77-0.54 


1.8(0.2) 


-3.5(0.8) 


9.3(1.4) 


93.1 GHz 








Orion KL 


3.5(0.3) 


7.6(0.3) 


3.2(0.3) 


NGC 6334F 


3.4(0.2) 


-8.3(0.4) 


3.6(0.4) 


94.5 GHz 








Orion KL 


14.7(0.4) 


7.6(0.1) 


4.5(0.1) 


NGC 6334F 


4.4(0.4) 


-7.9(0.3) 


5.0(0.6) 



emission is from a site north of the maser, but still within the 
telescope beam (Caswell et al. 2000). It seems likely that the 
85.5-, 86.6- and 86.9-GHz methanol transitions are exhibit- 
ing weak maser emission, while no corresponding thermal 
component was detected. 
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Radial Velocity (km s" 1 ) 

Figure 1. Methanol spectra observed at Mopra for Orion KL. 
3.3 NGC 6334F 

Methanol emission was detected at five frequencies in 
NGC 6334F (G351 .42+0.64), as shown in Fig. | The lines 
are broad, with no narrow features of distinctive maser 
character. An earlier detection of the 86.6- and 86.9-GHz 
methanol lines at low velocity resolution using the IRAM 
30-m telescope (Bachiller & Cernicharo, 1990) is consistent 
with our 3-Jy detections if the source is small with respect to 
the 25" IRAM beam, e.g. < 6" as found for the 156.6-GHz 
thermal line by Caswell et al. (2000). 

Maser lines in this source typically appear at velocity 
near -10 kms" 1 (Caswell et al. 1995a,b; 2000; Menten & 
Batrla 1989, Haschick et al. 1989), while thermal emission 
is slightly offset near -8 kms" 1 (Caswell et al. 2000, Slysh 
et al. 1999), as illustrated in Fig. ^ for the 6.6-GHz maser 
emission and the 85.5-GHz thermal emission. In our new 
observations all the methanol lines are near -8 kms -1 , and 
so all lines are probably of thermal origin. Because of the 
overlap between the thermal and maser velocity ranges, we 
use the peak intensity of the thermal line as an upper limit 
on the maser flux in the modelling which follows. 
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Figure 2. Methanol spectra observed at SEST (85.5 GHz) and 
Mopra (86.6, 86.9 GHz) for G345.01+1.79. For comparison we 
include spectra at 107.0 and 108.8 GHz from Val'tts et al. (1999), 
and at 6.6 and 12.1 GHz from unpublished observations made 
with the University of Tasmania 26m Mt Pleasant antenna. 



3.4 G351.77-0.54 

Weak broad methanol emission near -3 kms -1 was detected 
at 86.6 and 86.9 GHz in G351.77-0.54 (Fig. %, and is likely 
to be of thermal origin. The survey of Caswell et al. (2000) 
found similar broad emission at 156.6 and 107.0 GHz centred 
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Figure 3. Methanol spectra observed at SEST (85.5 GHz) and 
Mopra (86.6, 86.9, 93.1,94.5 GHz) for NGC6334F, with additional 
spectra at 107.0 GHz, 6.6 GHz, and 12.1 GHz, as in Fig. || 
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Figure 4. Methanol spectra observed at SEST (85.5 GHz) 
and Mopra (86.6, 86.9 GHz) for G351. 77-0.54, G9.62+0.20 and 
G10.47+0.03. Note the compressed velocity scale for G351.77- 
0.54. 



on -2.9 km s 1 in this source, which was the second strongest 
thermal source in their sample. 



3.5 G9. 62+0. 20 

The 85.5-GHz spectrum of G9. 62+0.20 shows a broad ther- 
mal component centred at 3.8 kms -1 and a narrow spike 
at -0.9 kms -1 , possibly of maser origin (Fig. |). Note that 
the spectrum of the 37.7-GHz methanol line (belonging to 
the same J_2 — (J + l)-i E line series) contains a maser 
component (Haschick et al. 1989) with velocity and width 
coinciding with that of the 85.5-GHz spike (to within 0.1 
kms -1 , which is less than the uncertainty 0.3 kms -1 in the 
rest frequency measurements) . Maser components with sim- 
ilar velocity are also present at 12.1- and 6.6-GHz (Caswell 
et al. 1995a,b). Although our 85.5-GHz spectrum is very 
noisy, the spike remains after substantial smoothing, and 
we regard it as a tentative detection of a maser component 
at the 3.6 sigma level. For comparison, the 107.0-GHz spec- 
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trum of this source displays two narrow spikes of maser emis- 
sion at -0.5 and 1.2 kms -1 with a broad underlying thermal 
component centred at 3.6 kms -1 (Caswell et al. 2000). 

3.6 G10.47+0.03 

Methanol emission at 3.3 Jy was detected at 85.5 GHz in 
G10.47+0.03 (W31), as shown in Fig. |. The broad line 
centred near 66 kms -1 is at the velocity typical of ther- 
mal emission for this source (see e.g. Caswell et al. 2000) 
and we conclude that the line is probably of thermal origin. 

4 MULTI-TRANSITION MASER MODELLING 
4.1 The maser model 

We employ the class II methanol maser model introduced 
by Sobolev & Deguchi (1994a) (hereafter the SD model). 
Methanol molecules are pumped by infrared radiation from 
warm dust to their second and first torsionally excited 
states, from which they decay, giving rise to population in- 
versions between certain levels of adjacent if— ladders in 
the torsional ground state. The extreme brightness tem- 
peratures displayed by some 6.6- and 12.1-GHz methanol 
masers are reproduced when the model includes amplifi- 
cation of free-free continuum radiation from an underlying 
ultra-compact Hn region, and when the spherical LVG ra- 
diative transfer approximation includes a 'beaming factor' 
which distinguishes the radial and tangential optical depths, 
to better represent an elongated maser geometry. Details 
of the pumping cycles are described in Sobolev & Deguchi 
(1994b). Model equations are summarized in the Appendix 
to Sobolev et al. (1997a), in which the SD model was first 
applied to A— species methanol. Maser predictions for many 
transitions under a range of representative model conditions 
appear in SCG97. 

The SD model contains three sources of continuum radi- 
ation: the infrared pumping radiation parameterized by Td, 
the cosmic microwave background, and the H n region free- 
free spectrum. The Hn region emits mostly at low frequen- 
cies and provides a source of background photons for am- 
plification by the low frequency masers. The Hn spectrum 
also influences the maser excitation, tending to counteract 
the IR maser pumping induced by the warm dust, via an in- 
crease in the number of low frequency photons penetrating 
the maser formation region. We vary only the geometrical 
dilution factor Whii, and use a fixed value for the turnover 
frequency of the Hn spectrum (12 GHz), and the electron 
temperature (the value of 18000 K used in the original SD 
model gives results very close to those obtained with the 
canonical value 10000 K). Radiative transfer is handled by 
the LVG approximation, with beaming factor e _1 defined as 
the ratio of optical depths parallel and perpendicular to the 
line of sight (Castor 1970). The tangential methanol column 
density Nm and the line width AV combine to give the pa- 
rameter Nm/AV, equivalent to the methanol volume density 
divided by velocity gradient. Maser brightness temperatures 
are largely determined by the optical depth along the line of 
sight, which is proportional to the product of Nm/AV and 
the beaming factor e _1 . 

The methanol excitation is also influenced by collisions 



with molecular hydrogen, parameterized by the gas tem- 
perature Tt and volume density riu- Unfortunately there 
are no accurately calculated state-to-state collisional excita- 
tion rates for methanol, and the empirical values we adopt 
(Peng & Whiteoak 1993) are based on a very small number 
of measurements. For comparison some calculations employ 
nonselective collisions, in which all downward rates between 
sublevels are assumed equal (scaled to give the same total 
cross- section) . 

Statistical equilibrium calculations in the SD model 
were carried out for a total of 1131 A— and E— species 
methanol energy levels up to J = 18 for v t = 0,1,2 
(Mekhtiev, Godfrey & Hougen 1999). Masers appear in 
many transitions over a wide range of values of each of 
the six model parameters described above, as sampled in 
SCG97. The parameter space is too vast to permit a full nu- 
merical exploration, so we identify trends by varying model 
parameters one or two at a time. Most results are based on 
an exploration of the cool gas regime with Tk = 30 K, typical 
of rotational temperatures derived from thermal methanol 
observations (Slysh et al. 1999). We adopt a dust tem- 
perature Td — 175 K which is sufficient to turn on the 
v t — 2 maser pumping, and within the range of tempera- 
tures identified for this type of source (Faison et al. 1998). 
Beaming factor e _1 = 10 and Hn region dilution factor 
Whii — 2x 10~ 3 are chosen to generate sufficient brightness 
temperature at 6.6 and 12.1 GHz to match VLBI observa- 
tions (> 10 K). The calculations span a range of parame- 
ters n H = 10 4 -10 9 cm" 3 and iViyr / A V = lO^-lO 13 cm~ 3 s. 
We have done additional less extensive calculations with 
e" 1 = 1, 3, 20, 100, Wmi = 10~ 2 - 10~ 5 , T d = 150 - 250 K, 
and Tk — 30 — 250 K. In general terms the presence or 
absence of individual masers in the model is governed by 
the values of the collisional parameters na and Tk, while 
their brightness is determined by the radiative parameters 
T d , Nm/AV, e- 1 and Whii. 

The class II methanol masers are all b— type transitions 
between states of adjacent K quantum number. Such transi- 
tions can be usefully grouped into series involving the same 
values of K but increasing or decreasing values of J. In 
the SD model many members of a transition series become 
masers under the same model conditions, consistent with the 
suggestion that some methanol K— ladders become overpop- 
ulated with respect to their neighbours (Wilson et al. 1985). 

The behaviour of some maser series is illustrated in 
Fig. ^ where we present the model dependence of maser 
brightness temperature on hydrogen density tih for selected 
methanol transitions in the cool gas regime. Many of the 
masers are present in the low density limit. Initial density 
independence is usually followed by growth of maser bright- 
ness, and eventual quenching as the density is progressively 
increased. No methanol masers are found at densities ex- 
ceeding 10 8 cm -3 . In the low density regime, which extends 
to approximately 10 5 cm -3 , the masers are generated by 
radiative processes and the collisions have very little influ- 
ence. The collisions are important in determining the peak 
maser brightness temperatures in the mid-density regime, 
particularly for those masers which are not present at low 
densities (e.g. 85.5 and 93.1 GHz). They also govern the 
point at which individual masers 'switch off' and hence the 
extent to which different masers coexist in the model. 

Some qualitative features of class II maser observations 
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Figure 5. Density dependence of selected class II methanol masers. Transitions are grouped into series and labelled with the line 
frequency in GHz, but for clarity not all members of each series are shown. The model parameters were as follows: kinetic temperature 
Tfe = 30 K, dust continuum temperature = 175 K, Nm/AV = 10 12,3 cm - 3 s, beaming parameter e~ 1 = 10, Hn region dilution factor 
Whii = 2 x 10~ 3 . 



are in accordance with Fig. kl The ubiquitous nature of 
the 6.6-GHz maser is consistent with its extreme bright- 
ness over a wide range of model conditions, while the less 
prevalent 107.0- and 156.6-GHz masers are weaker by or- 
ders of magnitude and have more stringent requirements 



for strong emission. In both modelling and observations the 
12.1-GHz maser is generally somewhat weaker than the 6.6- 
GHz maser. The same applies to lower J members of the cor- 
responding series, with the 108.8-GHz maser likewise weaker 
than its 107.0-GHz counterpart. 
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N M /AV=10 12 , e 1 =10, W H|| =2.10- 3 N M /AV=10 12 , e 1 =10, W H|| =3.10- 5 N M /AV=10 13 , e 1 =1 , W H|| =2.10 3 




HYDROGEN DENSITY / cm" 3 

Figure 6. Density dependence of the Vt = 0, Ji — ( J + l)o A + maser series for a model with e~ x = 10 (left panel), representing a beamed 
or elongated geometry. The centre panel shows the effects of a very distant H II region, while the right panel shows the effects of removing 
the beaming factor (spherical geometry). All calculations have Ty. = 30 K, = 175 K. 



Fig. |6j illustrates the influence of the Hn region back- 
ground continuum for two values of the geometric dilution 
parameter Whii- The H II region brightness temperature falls 
away sharply above 12 GHz, so it is not directly involved in 
the maser pumping, and its influence is mainly confined to 
the 6.6- and 12.1-GHz masers. These two transitions exhibit 
the strongest inversions of all the methanol maser lines, and 
are the first to undergo saturation. A reduction in VFhii (i.e. 
a more distant or smaller H II region) reduces the maser sat- 
uration by reducing the amount of low frequency radiation 
entering the maser region, while still providing a source of 
photons for amplification along the line of sight. The result 
is that the lowest frequency masers are greatly enhanced, 
while at millimetre wavelengths the effects are small. 

Fig. fj also shows the effects of changing the beaming 
factor e~ from the case where the maser region is opti- 
cally elongated towards the observer to the optically spher- 
ical case where the optical depth is independent of angle. 
The line-of-sight optical depth scaling factor e _ iViyr/AV 
has the same value for all models displayed in Fig. |(| and 
these models therefore have the same brightness in the ther- 
mal lines (i.e. those with excitation corresponding to LTE). 
For static clouds that would correspond to the same line-of- 
sight methanol column density. The beaming factor e _1 = 10 
represents an optically elongated maser region, in which the 
radial optical depth (along the line of sight) exceeds the 
tangential optical depth by a factor of 10. In this case the 
maser pumping is very effective, because much of the inci- 
dent infrared radiation encounters the lower tangential op- 
tical depth, allowing many pumping photons to penetrate 
the maser formation region. In comparison, in the optically 
spherical case with beaming factor e _1 = 1 (with a corre- 
sponding increase by a factor of 10 in Nm/&V to preserve 
the thermal line brightness), the optical depth encountered 
by the pumping photons is equally large in all directions. 



As shown in Fig. Q the resulting masers are much weaker, 
particularly at higher frequencies. This is attributed to the 
larger tangential optical depth in the spherical case, which 
limits the infrared pumping. 

4.2 Comparison of model with observations 

A difficulty in comparing models (with predicted bright- 
ness Tfc) with observations (with measured flux density) is in 
converting flux density measurements into maser brightness 
temperatures, since the required source angular size is too 
small to be measured with a single dish. We assume that the 
masers at different frequencies derive from a common area 
which is much smaller than the telescope beam, so their 
relative flux densities can be converted to relative bright- 
ness temperatures irrespective of the actual maser spot size. 
There is evidence that the class II methanol masers at dif- 
ferent frequencies are, at least, spatially coincident. The 
Menten et al. (1992) VLBI observations of the methanol 
masers in W3(OH) have absolute positional information for 
both the 6.6- and 12.1-GHz spots, and the strongest features 
coincide in both position and velocity. ATCA mapping of 
the 6.6- and 12.1-GHz masers by Norris et al. (1993) found 
an unambiguous coincidence of maser features within 0.02" 
in several sources. Low resolution interferometric observa- 
tions of W3(OH) at 23.1 GHz (Menten et al. 1988) and Cep 
A at 107.0 GHz (Mehringer, Zhou & Dickel 1997) demon- 
strate that these sources are small and have counterparts 
at 12.1 GHz. The 107.0-GHz masers were found to be co- 
incident with their 6.6-GHz counterparts to within typical 
measurement uncertainties of 5" for 8 sources checked by 
Caswell et al. (2000). 

Under the assumption that the masers are coexten- 
sive, the observed flux densities S(v) lead to the bright- 
ness temperature ratio Rots = [S(v)/ S(vq)\{vq/v) 2 with re- 



© 0000 RAS, MNRAS 000, 000-000 



Multi-transition study and new detections of class II methanol masers 9 



spect to a reference line at frequency vo. This may be com- 
pared with the model brightness temperature ratio R ca ic = 
Tb(v) /Tb(vo) . The advantage of comparing such ratios is 
that Tb{yo) may range over orders of magnitude between 
different models, in each of which a different source size 
is implicit. Brightness temperatures are constrained where 
VLB1 data is available. 

There are many reasons why a comparison between the 
model and a diverse collection of single dish observations 
can be at best semi-quantitative: the size of the maser re- 
gion is unknown, the masers at different frequencies may 
not be coincident, they may differ in extent, they may be 
variable, the observed positions are not always identical, 
there may be significant calibration discrepancies, etc. We 
estimate that values of R b s (which range over 5 orders of 
magnitude) may be uncertain by up to a factor of 3. Fur- 
thermore the current LVG model is likely to be deficient in 
reproducing details of the spectrum and distribution of the 
pumping dust, source structure, possible maser saturation 
effects, etc., and although radiative processes dominate the 
excitation, many quantitative details at high densities de- 
pend also on the poorly known collisional excitation rates. 
In view of such uncertainties, order-of-magnitude agreement 
between R b s and R ca ic may be all that can be achieved. 
Nevertheless, the maser brightness temperatures are so sen- 
sitive to the model parameters (Figs, g,^) that estimates of 
the hydrogen density and methanol column density can be 
obtained. 

4.3 G345.01 + 1.79 

Narrow lines at velocity near -21.8 kms -1 have been ob- 
served in G345.01+1.79 at 13 methanol frequencies, as sum- 
marised in Table ^. All these transitions become class II 
masers under some circumstances in the SD model (SCG97). 
This source is unusual in that the flux density of the 12.1- 
GHz maser at -21.8 kms -1 is not much less than that of 
the 6.6-GHz maser (Caswell et al. 1995b). Furthermore it 
has the brightest masers yet detected at 107.0, 156.6 GHz 
(Caswell et al. 2000) and the 157-GHz series (Slysh et al. 
1995), and is the only source in which narrow lines have 
been detected to date at 85.5, 86.6, 86.9 and 108.8 GHz. 
Upper limits for methanol emission at 93.1 and 94.5 GHz 
are also tabulated. 

In the following analysis we explore the ability of the 
SD methanol maser model to reproduce observational prop- 
erties of the G345.01+1.79 masers. We assume that all 13 
of the lines are of maser origin, and that they are formed 
in the same volume of gas. VLBI at 6.6 GHz has shown the 
maser spot size to be < 0.0036" with corresponding bright- 
ness temperature T 6 (6.6) > 10 11 ' 5 K (Phillips 1998), while at 
12.1 GHz the spot size is < 0.0007" and T 6 (12.1) > 10 110 K 
(Ellingsen et al. 1998). There is no mapping data available 
for the remaining transitions in G345.01+1.79, and evidence 
of small source size is required to fully establish their maser 
character, but similarities in maser line profiles suggest that 
the assumption of coincidence is reasonable. The masers in 
Table ^| are ranked by their brightness temperatures relative 
to that of the 6.6-GHz line. If the 6.6-GHz maser has bright- 
ness temperature > 10 1 K as suggested by the VLBI ob- 
servations, then under these assumptions the weakest masers 
detected have brightness > 10 7 K. 



We have compared the ratios Robs with those of the 
model calculations described in Section 4.1, with warm dust 
(Td — 175 K), cool gas (Tk = 30 K), and elongated ge- 
ometry (e" 1 = 10). The relative brightness of the low 
and high frequency masers in this source is approximately 
matched with a geometrically dilute H II region spectrum, 
Wmi = 2 x 10" 3 , giving #(107.0) ~ 10" 3 . General agree- 
ment for many lines is obtained for nu < 10 7 cm -3 and 
10 11 ' 5 < N M /AV < 10 12 - 5 cm _3 s. The model brightness 
temperatures at 6.6 and 12.1 GHz attain the VLBI values 
at the higher ends of these ranges. The difficulty of fitting 
a single component cool gas model to all the observations 
simultaneously is apparent from the trends in Fig. |5|. The 
86-GHz lines require relatively low densities, while the 85.5- 
GHz line requires moderately high density. There is some 
overlap between the 85.5- and 86-GHz masers at interme- 
diate density when the methanol column density is large. 
However the 94.5-GHz nondetection adds a further strong 
constraint, as it is difficult to find conditions in which the 
model predicts that the observed 85.5-GHz maser should be 
brighter than the undetected 94.5-GHz line. 

The results of some illustrative models are included in 
Table ^[ Model A represents the low density temperature- 
independent regime with Nm/AV = 10 12 ' 3 cm _3 s. It ac- 
counts for many of the observations, but is too bright at 
86 GHz, and predicts absorption rather than maser emis- 
sion for the 85.5-GHz line. Whether or not the masers are 
really coincident, the 85.5-GHz observation is clear evidence 
for a higher gas density origin for at least some of the masers. 
Model B is an intermediate density model (nrr = 10 6 cm -3 ) 
which matches the observed maser ratios reasonably well 
(including the 85.5-GHz line), though still too bright at 
86.9 GHz. However, it predicts maser action at 94.5 GHz 
comparable with that at 107.0 GHz, contrary to our ob- 
served upper limit which is more than an order of magni- 
tude fainter. Model C has lower column density of methanol 
(Nm/AV = 10 11 ' 8 cm _3 s) and gives a better prediction of 
the 86-GHz brightness while satisfying both the observed 
upper limits, but has the 85.5-GHz line in absorption. The 
difficulty of reconciling the 85.5-GHz maser observation with 
the 94.5 nondetection persists over the variety of model con- 
ditions which we have investigated, including higher gas ki- 
netic temperatures. Model D illustrates how a warmer gas 
temperature (Tk = 100 K), with a slightly higher hydrogen 
density and methanol column density, gives very similar re- 
sults to the cool gas Model C, again accounting well for all 
observations apart from the 85.5-GHz emission. 

When the calculations over a range of nn and Nm / A V 
are repeated with non-selective collision rates, the density 
at which some masers turn on or off alters by up to a factor 
of 30. For example, the 86-GHz masers extend to slightly 
higher density, while the 94.5-GHz maser is restricted to a 
lower range of densities than before. While the same range 
of riH and Nm/AV gives broad agreement with the obser- 
vations, the different collisional rates lead to changes in 
the details of conditions under which various masers co- 
exist. The best agreement with observations is exhibited 
by Model E of Table ^, with moderately large values of 
density (n H = 10 6 ' 5 cm~ 3 and N M /AV = 10 12 cm" 3 s), 
in which all brightness temperature ratios agree to within 
an order of magnitude and the upper limits are satisfied. 
Irrespective of the collision model, to obtain approximate 
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Table 4. Comparison of observations and modelling for class II methanol masers in G345.01+1.79. T;,(6.6) is the model brightness 
temperature of the 6.6-GHz maser, and model parameters are defined in Section 4.1. The relative brightness R oas is the ratio of observed 
(S/u 2 ) to its value for the 6.6-GHz line, while for the models we tabulate R ca i c = ' T b( !y )/Tf)(6.6) for comparison (see text Section 4.2). 
Model E uses non-selective collisions. 
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agreement between the maser model and the majority of 

the G345.01+1.79 maser observations requires parameters 

in the range n H = 10 5 ' 5 - 10 7 cm -3 and N M /AV = 
10 n.5 _ 1Q i2.5 cm -3 g with e -i = 10j while Tk ig not un _ 

ambiguously constrained. We conclude that more reliable 
collisional excitation rates are required to establish tighter 
limits on the model parameters. 

4.4 NGC 6334F 

Methanol masers at 8 frequencies have been identified in 
NGC 6334F near -10.4 kms" 1 , as summarised in Table ^| 
Identification of weak masers in this source is complicated 
by overlap with broad thermal emission, and we have taken 
the thermal flux density as a measure of the upper limit for 
the maser. This source was the second strongest 6.6-GHz 
maser in the 245 reported by Caswell et al. (1995a). VLBI at 
6.6 GHz has shown the maser spot size to be < 0.015" with 
corresponding brightness temperature T(,(6.6) > 10 10 ' 1 K, 
while at 12.1 GHz the spot size is < 0.0025" and r b (12.1) > 
10 10 ' 9 K (Ellingsen 1996). Imaging of the -10.4 kms" 1 emis- 
sion at 6.6 and 12.1 GHz shows that it is in the cluster 
coincident with the brightest continuum emission. 

Some of the differences between the methanol masers in 
NGC 6334F and G345.01+1.79 may be attributed to con- 
fusion by thermal emission, and absence of some key ob- 
servations, but the prevalence of high frequency masers in 
G345.01+1.79 and their general absence in NGC 6334F in- 
dicates significant differences in the pattern of maser excita- 



tion. The 6.6- and 12.1-GHz masers in NGC 6334F have sev- 
eral times greater flux density than those in G345. 01+1.79, 
while the millimetre wavelength masers are significantly 
weaker, such that only the 107.0-GHz maser has been de- 
tected in NGC 6334F, and this is weaker both in absolute 
and relative terms. NGC 6334F harbours 5 masers below 
40 GHz which have not been sought in G345.01+1.79. Three 
of these (37.7, 38.2, 38.4 GHz) belong to the same transition 
series as the millimetre masers seen only in G345.01+1.79 
(85.5, 86.6, 86.9 GHz respectively). 

In the following analysis we consider how the observa- 
tional properties of the NGC 6334F masers may be repro- 
duced within the SD model. The Hn region dilution factor 
Wmi is the model parameter with greatest influence on the 
maser brightness at low frequencies, via saturation effects 
in the strongest masers (see Fig. ^]and Section 4.1). Choos- 
ing Wmi = 3 x 10~ 5 reproduces the relative intensities of 
the 6.6- and 107.0-GHz masers, fl(107.0) ~ 10~ 5 . Apart 
from this smaller value of Whii, Model F in Table has 
parameters in the range found suitable for G345.01+1.79 
(n H = 10 60 cm" 3 , Nm/AV = 10 12 ' 2 cm - 3 s, e _1 =10). This 
accounts fairly well for all the observed masers except for 
that at 19.9 GHz, but fails to meet the upper limits estab- 
lished for the 86- and 94. 5- GHz lines. Model G has a smaller 
methanol column density, and satisfies the upper limits at 
millimetre wavelengths, while giving poorer agreement with 
the 20-40 GHz maser observations. 

An alternative to the very dilute Hn region spectrum, 
which is equally able to reproduce the trend of reduced 
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Table 5. Comparison of observations and modelling for class II methanol masers in NGC 6334F. Tf,(6.6) is the model brightness 
temperature of the 6.6-GHz maser, and model parameters are defined in Section 4.1. The relative brightness R ooa is the ratio of observed 
(S/u 2 ) to its value for the 6.6-GHz line, while for the models we tabulate R ca i c = Tf>M/T{j(6.6) for comparison (see text Section 4.2). 
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maser brightness at higher frequency, is to remove the beam- 
ing factor by setting = 1 (see Fig. |^). In this case the 
masers arise in a spherical region, with large column density 
of methanol, and hence optical depth, both laterally and 
along the line of sight (iV M /AV = 10 13 cm" 3 s). Model H 
illustrates the agreement with the observed maser bright- 
ness ratios, while Model I again shows that the millimetre 
wavelength upper limits are better satisfied with a smaller 
methanol column density. The reduced beaming factor re- 
sults in a smaller brightness temperature for the 6.6-GHz 
maser and a correspondingly larger source size, though still 
meeting the limits set by VLBI. 



5 DISCUSSION 

In comparing the model with a variety of maser observa- 
tions in two class II methanol maser sources we find that 
even in the best case the model predictions differ from the 
observations by as much as a factor of 10 for one or more 
transitions. We estimate that uncertainties in the measured 
flux densities may contribute up to a factor of 3 to this 
error, and the other major factor is likely to be inaccu- 
rate collisional excitation rates. A multi-component model, 
which relaxes our assumption that all the masers are spa- 
tially coincident, could potentially account for discrepan- 
cies in lines which are underpredicted, if their emission can 
be attributed to a region with different physical parameters 
(e.g. the 19.9-GHz maser in NGC 6334F, and the 85.5-GHz 
maser in G345.01+1.79). It will not, however, assist with 
the overprediction of maser emission at 94.5 GHz, accom- 
panying the maser detected at 85.5 GHz in G345.01+1.79, 



and this we attribute to the poor collision model. Despite 
the imperfect agreement between model and observations, 
order-of-magnitude estimates can be obtained for the hy- 
drogen density and methanol column density in the regions 
responsible for the majority of the maser lines. 

For both sources investigated, modelling suggests a 
gas density of approximately 10 6 cm -3 . Nm/AV is near 

10 12 cm _3 s with beaming factor e" 1 = 10, or near 

10 13 cm~ 3 s with e _1 = 1. For a typical maser line width 
AV ~ 0.5 kms -1 , this implies a radial methanol column 
density of 5 x 10 17 cm -2 in both cases. This is greater than 
values determined in the vicinity of class I methanol masers 
in the star formation regions SgrB2 and DR21 (Mehringer 
& Menten 1997, Liechti & Walmsley, 1997). Since the two 
sources studied are both exceptional examples of class II 
methanol maser sites (G345.01+1.79 being the strongest 
source of millimetre methanol masers, and NGC 6334F host- 
ing one of the strongest 6.6-GHz masers), it is reasonable 
that their methanol column densities represent the upper 
extreme. 

High resolution observations such as those by Caswell 
(1997) and by Phillips et al. (1998) show that the maximum 
extent of the maser spot distribution at a single maser site 
is approximately ~ 6000 AU (10 17 cm or 30 milliparsec). If 
we assume that we do not lie along a preferred line of sight, 
then the linear size of the maser cluster gives a reasonable 
estimate of the size of the masing gas cloud, and hence the 
path length. These assumptions suggest that the fractional 
abundance of methanol, Xm, in the masing regions exam- 
ined is ~ 5 x 10~ 6 . This is a little larger than the highest 
abundance observed by Mehringer & Menten (1997) toward 
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quasi-thermal cores in SgrB2 (4 x 1CT 6 ), but as the lat- 
ter is an average over a whole core, there should be regions 
of higher abundance. Such high methanol abundances are 
not predicted by existing chemical models. However, recent 
observations of a number of high mass protostars confirm 
earlier indications that methanol is the second most abun- 
dant solid-state molecule after water (see Dartois et al. 1999 
and references therein). This supports the hypothesis that 
the methanol abundance in the gas phase is enhanced by 
evaporation of methanol-ice from grain-mantles during the 
passage of a shock wave, and if this can be achieved without 
the simultaneous injection of gas then Xm ~5x 1CF 6 may 
be plausible. 

Class II methanol masers in the two sources modelled 
are consistent with maser regions exposed to radiation from 
warm dust (Td = 175 K). While the masers are generally pre- 
dicted to be brightest for cool gas (T fe = 30 K), higher gas 
temperatures are not inconsistent with the majority of the 
observations, and temperatures above the sublimation tem- 
perature of methanol ice (80 K) can more readily account 
for the high methanol abundances required. There are sev- 
eral possible reasons for the remaining discrepancies between 
model and observations, e.g. that the model has not been ad- 
equately explored in all parameter ranges, or that it is defi- 
cient in some important aspect of the treatment of radiative 
or collisional effects. In particular, since the collisional ex- 
citation rates are not well established, some aspects of the 
temperature or density dependence of individual masers are 
likely to be poorly modelled. Our simple model does not ac- 
count for the source structure, physical parameter gradients, 
or details of the dust spectrum. A further possibility is that 
our underlying assumption that all the masers are coincident 
is incorrect. High resolution observations will help to distin- 
guish these alternatives, as will the observation of further 
maser lines. For example the 80.9-GHz 72 — 8iA~ transition 
is predicted to be a significant maser in cool gas models (see 
SCG97), but not for gas temperatures Tk > 150 K. 

Since the observations were made at different times, 
some of the discrepancies between model and observations 
may be due to maser variability. Caswell, Vaile & Ellingsen 
(1995) found that 12.1-GHz masers were more variable than 
6.6-GHz masers in the same source, and that the weaker 
masers show greater variability. This was attributed to lower 
levels of saturation in the 12.1-GHz and weaker masers. We 
may therefore expect some variability in the weaker mil- 
limetre masers, and there is convincing evidence of this at 
107.0 GHz in one source (Caswell et al. 2000). The difficulty 
in accounting simultaneously for the 85.5- and 94.5-GHz ob- 
servations in G345.01+1.79 may perhaps be due to changes 
in the maser output over the intervening 2 years. The 20- 
40 GHz maser observations in NGC 6334F were made a 
decade earlier than the millimetre observations. 

The methanol maser spectra in the two sources dif- 
fer significantly, with the low frequency masers favoured in 
NGC 6334F while the high frequency masers are enhanced 
in G345.01+1.79. Modelling suggests at least two possible 
explanations: a more spherical geometry for the maser re- 
gion (beaming factor e _1 = 1), or a more distant Hn region 
(geometric dilution factor 3 x 10~ 5 ) in NGC 6334F. Since 
the model brightness temperatures of the masers in these 
two cases differ by two orders of magnitude, future high res- 
olution observations will help to elucidate the situation. 



6 CONCLUSION 

We have made the first detections of weak class II methanol 
masers at 85.5, 86.6 and 86.9 GHz; these are predicted by 
a pumping model involving torsional excitation to v t = 2 
by infrared radiation from warm dust. By assuming that 
the masers at different frequencies are coincident we have 
made the first attempt at multi-transition maser modelling 
for methanol. High resolution observations at many frequen- 
cies are required to explore these assumptions. Further ob- 
servations of the established masers in the frequency range 
20-40 GHz would also be valuable for future studies of this 
type. 

The level population and radiation transport mod- 
elling of G345.01+1.79 and NGC 6334F suggests that their 
methanol masers form under similar conditions of hydro- 
gen density near 10 6 cm~ 3 and column density of methanol 
as large as 5 x 10 17 cm~ 2 . The physical parameters implied 
from the maser modelling of these sources are generally con- 
sistent with those determined from observations of massive 
star formation regions. The major differences between the 
maser spectra in the two sources can be attributed to differ- 
ences in the geometry of the maser region, e.g. beamed or 
elongated geometry favouring the development of the high 
frequency masers in G345.01+1.79, while the dominance of 
the lower frequency masers in NGC 6334F is consistent with 
a more spherical maser region. 

The large number of different methanol maser lines ob- 
served in these two sources allows us to make a comparison 
with the simple SD maser pumping model, and to identify 
some particular shortcomings. Among these is the model 
prediction of significant maser action at 94.5 GHz, exceeding 
our observed upper limit, accompanying the newly detected 
85.5-GHz maser. We conclude that the 94.5-GHz maser is 
quenched in the maser sources, but that our current best es- 
timates of collision cross-sections do not adequately describe 
the excitation of this transition. The density and tempera- 
ture determination is imprecise because of such uncertain- 
ties. 

We have identified the ratio of high and low frequency 
maser intensities as a key feature distinguishing the two 
sources studied in detail. This approach is likely to be use- 
ful in the analysis of other methanol maser sources, since 
there is evidence that the 157-GHz maser series has been 
detected most readily in sources which are also strong in the 
107.0-GHz maser. Of the 13 sources which harbour masers 
at 107.0 GHz and were searched at 157 GHz (Val'tts et al. 
1995, Caswell et al. 2000, Slysh et al. 1999), the four sources 
in which 157-GHz masers were identified (Slysh et al. 1995) 
are among five with 107.0-GHz masers exceeding 20 Jy. 

A multi-transition study of methanol masers requires 
many assumptions, and even the relatively simple model 
adopted here has a vast parameter space which has been only 
partially explored. Nevertheless, it is encouraging that semi- 
quantitative trends can be identified from the data presently 
available, and an explanation of the intensities of all the 
class II methanol masers seems to be possible in principle. 
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